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Abstract 

Maximal helicity- violating scattering amplitudes in = 4 supersymmetric Yang-Mills theory 
are dual to Wilson loops on closed null polygons. We perform their operator product expansion 
p I analysis in two-dimensional kinematics in the soft-coUinear approximation which corresponds to 
D . the case when some light-cone distances vanish. We construct the expansion in terms of multi- 



particle "heavy" -light operators, where the "heavy" fields are identified with the Wilson lines 
defining the OPE channel and the light fields emerge from the curvature of the contour. The 
correlation function of these define the remainder function. We study the dilatation operator 
for these operators at one loop order and find that it corresponds to a non-compact open spin 
chain. This provides an alternative view on elementary excitations propagating on the GKP 
string at weak coupling, which now correspond to particles traveling along an open spin chain. 
The factorized structure of the Wilson loop in the soft limit allows one to represent the two- 
loop correction to the octagon Wilson loop as a convolution formula and find the corresponding 
remainder function. 



Maximal helicity-violating (MHV) scattering amplitudes in maximally supersymmetric Yang- 
Mill theory are dual to Wilson loops stretched on null polygonal contours. The latter is built 
from the on-shell momenta of massless scattered particles. The equivalence of the two objects was 
hinted at strong coupling fl6\ and then established at weak coupling at one |2l [3] and two-loop 
orders with up to six sides [2J putting the conjecture on a firmer foundation. Dual conformal Ward 
identities |4] allow one to constrain the form of the Wilson loops with an arbitrary number of sides. 
These arise from the known pattern of the dual conformal symmmetry breaking due to ultraviolet 
divergences of loop integrals from the presence of cusps on polygonal contours reproducing the 
ABDK/BDS ansatz A conformally invariant remainder function, which can only depend on 
conformal cross-ratios of light-cone distances, enters as an additive component of the full answer 
for the Wilson loop starting from the two-loop order and involving more than five segments on its 
contour. It is the main unknown in confronting both sides of the duality. Its computation in the 
simplest case of the hexagon serves as a very illustrative example which exhibits inefficiency of the 
ordinary techniques based on Feynman integrals. First, agreement was established numerically 
between the corresponding six-leg amplitude [6] and Wilson loop |7] expressions. The follow-up 
tedious analyses culminated in analytical results for the latter filling dozens of pages in terms of 
Goncharov polylogarithms [S]. Eventually, the theory of motives [H] allowed one to reduce the 
entire output to a two-line expression in terms of classical polylogarithms [TU]. Hence, as this 
example clearly demonstrates, explicit perturbative computations are extremely complicated, 
hindering the unraveling of a simple internal structure of observables and thus techniques that 
can bypass involved loop analyses, getting to the final answer with an efficient shortcut would 
be highly valuable. A little while ago, a formalism based on the operator product expansion was 
suggested in Ref. [H] and used to compute the two-loop remained function for polygons in the 
restricted two-dimensional kinematics [I2] as well as the hexagon Wilson loop in full R^'^ [I3] 
in a very concise manner. Quite recently, this technique was also used to constrain the symbol 
[21 [To] of the three-loop hexagon [H] . 

In this note, we provide a complementary view on the operator product expansion by ex- 
plicitly constructing the first contribution to the series in terms of "heavy" -light operators and 
demonstrating how their renormalization group evolution determines the expectation value of the 
Wilson loop. This consideration also provides an interpretation for Gubser-Klebanov-Polyakov 
(GKP) excitations [151 [HI [IZ] from the point of view of open spin chain^. We limit our consid- 
eration to the two-loop remainder function and thus will not (need to) address loop corrections 
to the coefficient functions of contributing light-cone operators, though this should not represent 
a problem of principle for the current systematic expansion. 

To set up the formalism, we restrict our analysis to the two-dimensional kinematics and 
thus address the simplest nontrivial case corresponding to the octagonal null contour. The object 
at the center of our analysis is the Wilson loop 

ly = ^(tr(iy[o,i]...%,8]W^[8,l])), (1) 

built from eight light-like segments 

Wiii+i] = T exp {ig dU x{ti) ■ A{U)^ , (2) 

^Unrelated open spin chains had made their appearance in the study of Regge hmit of scattering amphtudes 
in Ref. [18]. 
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running along particle's momenta as x{ti) = Xi — tiXu+i. Below, we will dress the Wilson line 
segments with the superscripts ± to designate the light-cone direction of the corresponding links 
on the zig-zag contour (as shown in Fig. [T]). The Wilson loop expectation value admits the 
perturbative expansion 

W = J2 , (3) 

n=0 

in the 't Hooft coupling a = g'^^Nc/{87!'^). Since we will be after the remainder function Wr 
depending only on the conformal cross-ratios, in the current case one can construct just two of 
these, which we conveniently choose as 

^+^V67^ = (4) 

Xg2X73 •^84"''73 

To start with, we will choose a channel for the operator product expansion as in Ref. [H], say 
in the x+ direction. We introduce Wilson lines for segments running long the top and bottom 
side of the loop as follows, Wtop = ^[56] ^[6 7] ^[7 8] ^hot = ^[T2]^[2 3]^[3 4]' respectively. 
Then, in the short light-cone distance limit x'^2 ~^ ^^"^ ^67 ~^ 0' when the segments adjacent 
to these become coUinear, the operator product expansion for these two sides can be constructed 
systematically [19J, with the leading terms having the form 

W^topL+^O = 47^[5,0]W^[0,8] + O {{Xt,f) , H^botL3+^o = ^2 ^^^^^0' ,1] + ^ ((^S^)') , (5) 

and being induced by the gluon field strength insertion into the contour, 
J'[,,o]='9 [ dtW^l^^x-{t)F^'{t)W^l^^, J'yo',i]=^9l dt'W^,,,^x'-{t)F^-{t')W^,^,y (6) 

•/ •-'0 

Here x{t) = xq — txoe and x'{t') = xq' —t'xo'2- Notice that the soft-collinear expansion eliminates 
several cusps of the original Wilson loop and thus becomes insensitive to the corresponding 
ultraviolet divergences. However, presently we are interested only in the remainder function which 
is free from the latter. Higher order terms in the series ([5]) can be constructed systematically 
in an analogous fashion and are expressed in terms of multi-field insertions along corresponding 
contours |T9]. 

The Wilson loop is determined by the correlation function of the above operators 

where the ellipses stands for higher order terms in x^2 ^'^d x^j. At leading order in 't Hooft 
coupling, this correlation function is simply determined by the one-gluon exchange (see the right 
panel in Fig. [1] where one sets Sj = 1) and immediately yields 

= -if^^o'Y^^^^ [i^m'T - i^oo'T + i^mT - #1 ^ log(l + U-) . (8) 

fc {Xqq, ) 

Here the last factor in the dimensionally-regularized expression was replaced (in the final equality) 
by the invariant cross-ratio to the accuracy that we work with since the result is valid for 
— )■ 0. In order to restore the full result, we choose another channel for the short light- 
cone distance expansion of the Wilson loop, i.e, involving the ^"-direction this time. Then, 
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Figure 1: Null Wilson loop in the R^'^ kinematics (left) and leading term contribution to its 
operator expansion in the soft-coUinear limit x^g, Xqj — )■ 0. 

performing the same analysis as done above, one finds analogously that W^^ = —u^ log(l + W^) 
as xj^4,Xg7 0. Bootstrapping the result, one immediately obtains the exact expression = 
— log(l + M") log(l + u~^) valid for arbitrary values of in agreement with |12] . 

Beyond the one-gluon exchange, one has to compute corrections to the above correlation 
function of the light-cone operators J-"". The calculation is split into the analysis of the one- 
loop corrections to coefficient functions and evolution kernels. By choosing the light-cone gauge 
A'^ = 0, one immediately observes that all gauge links on the top and bottom of the loop 
degenerate to unity, i.e., W~ = 1. Therefore, the only interaction left is of the exchanged 
F"'"~-field with the Wilson loop segments running long the direction. 

Since we can use the OPE in different channels to reconstruct the full function by means of 
the bootstrap procedure, we can constrain the -u'-dependence at two loop order by studying the 
large evolution x[|Q,-time limit of the correlation function. The leading effect will arise from the 
logarithmically enhanced terms. As we can already see from the leading order result the time- 
evolution is associated with the dimensionless combination /ixjg,. Notice that the dependence on 
the factorization scale /i will cancel upon accounting for the short distance coefficient function 
which will introduce logarithms involving short light-cone distances entering in dimensionless 
combinations fix^2 bottom portion of the loop and fJ^x^j for the top one, thus restoring 

the complete cross-ratio in the sum. This is completely analogous to the conventional DGLAP 
logarithms in the analysis of structure functions of deep-inelastic scattering. Thus, all we have 
to do is to calculate the leading order dilatation operator for the effective light-cone operator 

<^[xi,X2,X3] = h^{x:[)F+-{x2 )h{x^) , (9) 

built from the gluon field strength and two "infinitely-heavy quark" field^ h(xi) = W^[^o] 
and h^xs) = W^^-^, with the paths starting and ending at the points Xq = (xf, L"*"), x^ = (xj , 0) 

^The analogy is drawn with heavy-quark effective theory where the heavy quark gets replaced by a Wilson 
line along its classical path and interacting with external gauge fields through eikonal vertices. The difference 
however is that heavy quarks trace time-like trajectories while we are dealing with light-like ones. However, 
due to superficial similarity of the two object, we will call the light-like Wilson loop arising in our consideration 
"heavy" since they will not recoil against massless fields building up light-cone operators similarly to their massive 
counterparts. 
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and xi = {x^ , 0), X4 = (X3 , respectively, and parametrized by the corresponding trajectories 
x{t) = xq — txoi and x(t) = X3 — txu- We set the length of the path L"*" equal to the +-length of 
the Wilson loop Xqq,. This will serve as a regularization parameter of the well-known light-cone 
singularity of the null Wilson line. We did not display the gauge links on the intervals [1, 2] and 
[2, 3] since they reduce to unity in the A"*" = gauge. 

The dilatation operator can be easily computed from the Feynman graphs shown in Fig. |2l 
Making use of the well-known coordinate-space perturbative techniques (reviewed, for instance, 
in Appendix G of Ref. [20]), one immediately finds the following result for the pole-part of the 
one-loop correction to $, 

^xi,X2,X3\ = (^1 + ^i3i23) '^'(o)[a;i,a;2,a;3] (10) 

where kernel of the dilatation operator has a pair-wise structure at this order 

iOl23 = ^Ii2 + i023 (11) 
with individual Hamiltonians taking the form 



(12) 



Notice the presence of the double-pole in e which is a familiar consequence of the ultra-violet 
divergences generated by integrations in the vicinity of the cusp points of two adjacent Wilson 
lineqj that build up the operator $. These contributions are accompanied by the one- loop cusp 
anomalous dimension. They will not be of interest for our subsequent consideration since they 
are exactly captured by the dual conformal Ward identities However, the second term in the 
above expression, which is the focus of our analysis, induces the displacement of the light field 
F~^~ along the light cone and encodes its renormalization group evolution. Intuitively, only the 
light degrees of freedom get shifted since the "heavy" fields encoded by the Wilson lines do not 
recoil against the former. The one- loop graphs displayed in Fig. [2] give 



[)i2$[xi,a;2,a;3] = / da—- — axi 0x2, X3] , (13) 
Jo fJ+ 

()23'^'[a;i,a;2,X3] = / da ^Xi,ax2 + 0x3^x3] , (14) 

Jo m+ 

with a = 1 — a. Here the calculation was performed for an arbitrary propagating field, — 
anticipating follow-up applications to general R^'^ kinematics and supersymmetric Wilson loop, 
— which is encoded in the value s of the conformal spin of the excitatioE0. For the case at hanql. 



•^One of the lines gets degenerated to unity so that the corresponding divergence becomes end-point. 

""^These are gluon field strengths classified with respect to their twist, for the bosonic Wilson loop, or scalars 
and fermions, for the contour lifted to superspace. 

^Notice that the twist-three excitations like are not quasipartonic. As a consequence, the explicit calcula- 
tion in the light-cone gauge v4+ = of the d^A^ — > d^A^ transitions has to be supplemented by d± A± 
ones, with the latter traded for the d~^A~ contribution by mean of equations of motion. Notice that we also ig- 
nored two-to-three transitions which do not contribute at this order in coupling to the remainder function in the 
restricted two-dimensional kinematics. 
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it is with s = 1. In the above dilatation operators we introduced a conventional notation 
for the regularized end-point singularity 

1 1 fU{3 
— o[a) 



[a]+ a Jq (3 

The pairwise Hamiltonian f) is not s/2 invariant, and commutes with the generators of the light- 
cone algebra up to boundary terms (for J~) 

The above three-site Hamiltotian f)i23 = f)i2 + f)23 is also invariant under three-site generators 
J"^ = J2i Jii+i to the boundary terms, 

[f)l23,J°]=0 [t)l23,J+]=0 [t)i23,Jl=i:-, (16) 

where L~ = xj~3 is the length of the operator $. 

Making use of this result, we can instantly extract the logarithmically enhanced correction to 
the Wilson loop at two loops. Being after the remainder function, we can then promote logx^Jg, 
to — llogM"*", which is the only conformal-invariant cross ratio involving the "time" -evolution 
interval. The contribution takes the form of pair-wise kernels of the dilatation operator acting 
on the "heavy" -light system that defines the OPE limit of the null Wilson loop, 

^L+,x+^0 = |4742log(/^4o') (17) 
X (tr (f)[W^[0',5]-^[5,0]]^[0,l]-^[T,0'] + ^[0',5]-^[5,0]^)[^[0,l]-^[T,0']] + (^0' ^ xo))) + . . . . 

Here the two-particle operator i) acts on the "heavy" fields W and the field strength tensor F~^~ 
in J-'~ according to Eqs. (fT3l) . In the above equation, we restored the broken conformal symmetry 
of the evolution kernels by adding a parity-reflected contribution xq' ^ xq which restores exact 
commutation relations of the s/2 algebra. The calculation is straightforward and can be simplified 
making use of the conformal symmetry of the expected result. Making use of the explicit action 
of the Hamiltonians on pair-wise systems of fields and choosing the external points as follows 
Xq, = 0, X2 = e~°", Xg = 1,Xq = 00 such that cr = log(l + l/M~), we obtain the following integral 
representation of the two-loop remainder function in the soft-collinear limit of the a;"'"-channel, 
which yields the final answer. 



H^f = -VlogM+ /" y^{log(l-ae~'^) + log(l-a-V'^)} (18) 

^ Jo N+ 

= — -n^ logM^ logu^ log(l + u^) . (19) 

Bootstrapping this result by considering the OPE in the cross x~-channel, we reconstruct the 
two-loop remainder function in the form wj^^ = —1 logu"*" log(l + u~^) logw" log(l + u~). This is 
indeed the expression obtained by explicit diagrammatic calculations in Ref. [21] and confirmed 
by the OPE techniques in Ref. [TT] (see also [22]). 

The above reformulation of the OPE for the Wilson loop in terms of "heavy" -light light-cone 
operators provides a new perspective on the relation of the Wilson loop to integrable structures. 
To make a connection to the OPE formalism elaborated in Ref. [TT], it is instructive to point out 
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Figure 2: One-loop dilatation operator for the "heavy" -light operator $. 

a relation of the current problem to non-compact open spin chain systems, see Refs. [231, 125]. 
Since the Hilbert space of the pair- wise Hamiltonian does not possess lowest/highest weight 
vector^, the Algebraic Bethe Ansatz techniques are not applicable to the problems at hand. 
However, the solution can be found making use of the Baxter 0.{u) operator. Namely, the two- 
particle Hamiltonian f) in Eg. ( 1T3|) can be deduced from which obeys the Baxter equation 
for a two-site open spin chaiiiLj 

(M-Zs)n(M-z) = t2(u)n(M), (20) 

where the local integral of motion enters through the function i2{u) = u — p with p = = 
i{xdx + s), which is the 22-component of the Lax operator Sl{u) = u + icr ■ J . Since the equation 
( 120|) is merely a one-term recursion relation, it can be easily solved with the result 

Qt") = A'M^^^. (21) 

1 [lU + s) 

where N{u) is an-i periodic factor, N{u — i) = N{u), which we choose to be one. Then the 
Hamiltonian is constructed from its logarithmic derivative as follows 



f) = z — 

du 



log - i){2s) + ^(1) = -^(s - ip) + ^/'(l) 

da 







a 



a-'^+'-^ - 1] , (22) 



where we added the anomalous dimension per excitation. This coincides with the hamiltonian 
computed at leading order of perturbation theory (fT3|) . Here the momentum operator p possesses 
continuous eigenvalues and encodes the momentum of the light field propagating along the chain. 

The three-particle Hamiltonian acting on the operator $ in Eq. (Q, which is defined on a 
semi-infinite interval, 



[)i23^'[0,X2,oo] = / -—^a<l>[0, 0x2,00] + / 7-— $[0,ax2,oo] , (23) 
Jo N+ Ji N+ 

can be easily diagonalized via the Mellin transform 

<l>[0,X2,oo] = / dpx-'P-''^, (24) 
c 



^This is analogous to the absence of local operator expansion for heavy-light wave functions of the heavy-quark 
effective theory [37]. 

''However, only one site is dynamical, i.e, the one corresponding to the light field. 
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such that the total Hamiltonian 



f)l23$[0,X2,Oo] 




dpx2 



ip—1 



e{jp) $, 



(25) 



possesses the eigenspectrum 



e{p) = V^(l - ip) + ^{1 + ip) - 27/^(1) , 



(26) 



parametrized by the momentum eigenvalue p. This reproduces the energy of a single hole propa- 
gating on the background of derivatives which describes the elementary GKP excitation at weak 
coupling [261 IH]- Making use of the operator product expansion for the null Wilson loop, we 
can diagonalize the two-loop correction to the remainder function ( [T7|) by performing the Fourier 
transformation in a (or Mellin transformation in X2 as above) of the leading order coefficient 
function (jH]) and obtain the representation 



which coincides with the one advocated within the OPE framework of Ref. |TT] for the null 
polygonal Wilson loop. Here, the contour C runs slightly above the real axis. In this equation, 
we reconstructed the ^'''-dependence exactly by bootstrapping this expression making use of the 
cross channel as was done above. 

The formalism proposed in this note provides an explicit framework for the OPE analysis of 
the null polygonal Wilson loops. It allows one to clearly separate loop effects stemming from 
coefficient functions and evolution kernels and thus suggests a systematic extension for addressing 
higher order perturbative corrections. Our consideration also proposes a novel interpretation of 
the GKP excitations from the point of view of a non-compact open spin chain rather than a closed 
Heisenberg magnet picture involved in the analysis of the large-spin asymptotics of anomalous 
dimensions of twist operators. In the latter, the GKP modes emerge as holes propagating on the 
background of derivatives building up Wilson operators. 

Recent applications of the operator product technique to the supersymmetric case [28] do not 
heavily rely on the pertinent super Wilson loop [29l |30] and thus do not directly elucidate its role 
in the relation to non-MHV amplitudes which currently is not defined beyond tree approximation 
due to lack of a supersymmetry preserving regularization scheme required for its robust analysis 
[31] . The method employed in Ref. [28] is not particularly sensitive to the specifics of the coupling 
of GKP excitations to the super-contour of the loop. We hope that the technique we advocated 
in this paper will allow one to probe the connection between the scattering amplitudes and 
super Wilson loops beyond MHV approximation in a straightforward fashion. It would also be 
interesting to study higher-loop corrections to the open-spin chain hamiltonian and establish a 
relation to the framework of factorized scattering of hole excitations on the large-spin closed 
noncompact spin chains. 

We would like to thank Benjamin Basso for interest in the project at its very early stage 
and Gregory Korchemsky for reading the manuscript and instructive remarks. This work was 
supported by the National Science Foundation under grants PHY-0757394 and PHY-1068286. 
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